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Resonance Fluorescence Squeezing Spectra from Photonic Bandgap Crystals 
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The fluorescence intensity and quadrature spectra from a two-level atom embedded in a photonic 
bandgap crystal and resonantly driven by a classical pump light are calculated. The non-Markovian 
nature of the problem caused by the non-uniform distribution of the photonic density of states is 
handled by linearizing the generalized optical Bloch equations with the Liouville operator expansion. 
Unlike the case in free space, we find that the bandgap effects will not only modify the fluorescence 
spectral shape but also cause squeezing in the in-phase quadrature spectra. 
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The study of fluorescence spectra from two-level atoms 
have been a central topic in quantum optics since the 
beginning era of quantum mechanics in 1930's. From 
the view point of light scattering, both elastic Rayleigh 
scattering and inelastic Raman scattering processes are 
involved [l| and thus the fluorescence spectra will have 
a triplet shape as first calculated by Mollow 0. Theo- 
retical calculations of the fluorescence spectra has been 
explored 0,0] and also verified in experiments 5]. The 
squeezing phenomena in the phase-dependent fluores- 
cence spectra of the quadrature field component were first 
predicted in 1981 by Walls and Zoller [6| and also in 1982 
by Mandel pfl. It has been theoretically shown that the 
squeezing can be found in the out-of-phase quadrature 
component spectra under the condition that il 2 < -|r 2 
where SI is the Rabi frequency and V is the atomic 
decay rate. Squeezed fluorescence spectra have also been 
experimentally observed in an experiment using 174 Yb 
atoms [To|. 

In the theoretical studies of fluorescence spectra for the 
free space case, the Markovian approximation is usually 
used to describe the statistical properties of the optical 
noises. This is a good assumption for the free space case 
but is not a ppl icable for the case of photonic bandgap 
crystals 0, . This is because in a photonic bandgap 
crystal, the propagation of light is prohibited within a 
certain range of wavelengths (the bandgap) due to the 
lack of available photon states. Because of the bandgap 
effects, the distributions of the photonic density of states 
(DOS) are typically highly non-uniform near the band- 
edge. Such a property has prohibited the direct appli- 
cability of the Markovian approximation to simplify the 
derivation for the problem we are going to consider. 

In recent years the atom-photon interaction in pho- 
tonic bandgap materials has been found to exhibit many 
interesting new phenomena 0, Q, ^| . But almost all 
of these studies are focused on the transient behavior 
of the atom-photon interaction and to the best of our 
knowledge there is still no theoretical treatment on cal- 



culating the steady state fluorescence squeezing spectra 
from photonic bandgap crystals. The aim of this paper is 
thus to investigate the properties of the steady state reso- 
nance fluorescence emitted by a two-level atom embedded 
in a photonic bandgap crystal and driven by a classical 
pumping light. We treat the photon states in the pho- 
tonic crystal as the background reservoir and obtain a set 
of generalized Bloch equations for the atomic operators 
by eliminating the reservoir field operators. The non- 
uniform DOS distributions near the bandedge are mod- 
eled by the three-dimensional anisotropic dispersion rela- 
tion [lq and the Liouville operator expansion is used to 
reduce the two-time atomic operator products into equal- 
time atomic operator products. In this way the nonlinear 
generalized Bloch equations are reduced into a set of lin- 
ear equations with memory function terms caused by the 
atom-reservoir interaction. This set of linear equations 
can then be easily solved in the frequency domain and the 
resonance fluorescence spectra can be directly obtained 
from the correlation functions of the atomic operators 
in the frequency domain. After performing numerical 
calculation, we find that the spectral shape of the fluo- 
rescence intensity spectra will vary with the wavelength 
offset between the atomic transition wavelength and the 
bandedge. More interestingly, squeezing phenomena are 
found to present in the in-phase quadrature spectra in- 
stead of the out-of-phase quadrature spectra. 

To begin the derivation, the Hamiltonian for the sys- 
tem to be considered can be written as: 



H = -m a a z + h^iVkal 



(1) 
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Here the transition frequency of the atom and the fre- 
quency of the pumping light are denoted by u> a and 
tor respectively, a\. and are the creation and anni- 
hilation operators of the photon states in the photonic 
bandgap crystals, SI is the Rabi-flopping frequency of 
the atom under the external pumping light and also 
represents the relative magnitude of the pumping light, 



(|2)(2|-|1)(1|), 
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ert are the usual 
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Pauli matrices for the two-level atom, and g k is the atom- 
field coupling constant. Here we have used the index k 
to label different photon states and the coupling constant 
g k can be expressed as: 



9k(d,f ) = g k = \d\ui a 



1 



2he Q u k V 



d-E^(r ) (2) 



Here we use notations |d| for the magnitude of the atomic 
dipole moment, d for the unit vector along the direction 
of the dipole moment, V for the volume of the quantiza- 
tion space, and eo for the Coulomb constant. 

Starting from the Hamiltonian, we treat the photon 
field operators as the background reservoir and eliminate 
their corresponding equations to obtain the following set 
of generalized Bloch equations. 



*-(*) 



a+{t) = -i-a z (t) 



iAI 



+ / dt'G(t-t')a z {t)a-(t')+n-{t) 



iAI. 



(3) 



(4) 



dt'G c (t - t')a + (t')a z (t) + n+(t) 
& z (t) = in(a-(t)e iAt - a + (t)e- iAt + n z (t) (5) 
2 ! dt'[G(t-t')<7 + {t)a-(t') + G c (t-t')<7 + {t')a-(t)} 



Here A = lul — w Q and Afe = u> a — uj k . The two 
functions G(t) and G c (t) are the memory functions 
of the system caused by the atom-reservoir interaction 
and they are defined as G(t) = J2 k \9k\ 2 e lAkT Q{T), and 
Gc(t) = J2 k |fffc| 2 e" lAfcT e(r). Here 9(r) is the Heaviside 
step function. The three noise operators rz_(t), n + (t), 
and n z (t) originate from the original photon filed opera- 
tor before interaction and are expressed as follows: 



n-(t) = i^2g k e iAkt a z (t)a k (-oo) 



(6) 



n + (t) = -iJ29* k e- iA * t al(-™)a z (t) (7) 

k 



n z {t) = -2iY,9ke lAkt <?+(t)a k (-oo) 



(8) 



-2i 



k 



Supposing that the reservoir is in thermal equilibrium, it 
can be easily shown that at zero temperature the three 
noise operators are zero mean with only four non-zero 
correlation functions, (n^(t)n+{t')) r, (n-(t)n z (t')} r, 
(n z (t)n + (t')) R , and (n z (t)n z (t')) R . 

To actually evaluate the memory functions as well as 
the correlation functions of the noise operators, one needs 
to know the spectral distribution of the photonic den- 
sity of states. Although in general the DOS of photonic 
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FIG. 1: Amplitude and phase spectra of the memory func- 
tion G(lo) with bandedge frequency u c = 100/3. The memory 
function is non-uniform around the bandedge and becomes 
pure imaginary inside the bandgap. 



bandgap crystals is very complicated and also varies with 
the geometrical structure and the dielectric constants of 
the material, it is still possible to approximately model 
the DOS near the bandedge with a simple formula. Ac- 
cording to the results from the full vectorial numeri- 
cal calculation, the DOS near the bandgap for three- 
dimensional photonic crystals increases from zero and 
behaves more like the anisotropic model proposed in the 
literature |l6[ . To be more specific, for a three dimen- 
sional photonic bandgap crystal, if the wavevector that 
corresponds to the bandedge is kg, then the dispersion 
relation in the anisotropic model is described by the fol- 
lowing form: Lo k = oj c + A\k — k | 2 , where A is a model 
dependent constant and lu c is the bandedge frequency. 



Based on this dispersion relation, the corresponding DOS 
is given by: D(ui) = A l/ 2 y/u> — oj c Q(uj — u c ). The mem- 
ory functions under this anisotropic model also can be 
derived as: 



UJ C + y/uj c - LO a - LO 



3/2_ 



(J c + y/cj c - UJa+U) 



(9) 
(10) 



r], and we have used the space av- 



where /3 3 / 2 

erage coupling strength r) = ^ J dfl\d-F,\ 2 m the deriva- 
tion. 

From Fig. ^ the spectrum of the memory function 
G{u>) is non- uniform and asymmetric as we expect. When 
the frequency is below the bandedge frequency uj c , the 
memory function becomes pure imaginary, indicating 
the inhibition of the spontaneous emission inside the 
bandgap. The spectrum for another memory function 
G c (lu) is also similar. It can be easily checked that 
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the full-width-half-maximum (FWHM) bandwidth of the 
memory functions in Eq. © and Eq. i|l(J|) are 4uj c . For 
the bandgap in optical domain, the order of uj c is about 
2.014-15 Hz, and the typical lifetime of the atom is from 
1CP 3 sec to 1CP 9 sec, which is much longer than the re- 
sponse time of the memory functions. 

The generalized Bloch equations are a set of non- 
linear operator equations and cannot be solved eas- 
ily. To overcome this difficulty, we introduce the 
following Liouville operator expansion: cry(i) = 



n= 0.25 p 
<o s = 1 oo p 



l-i(t-t')V 



■L n a lj (t'). Here the 



Liouville super operator C is defined as C n <Jij{t') — 
7pr[- • • [<7ij(t'), H],H], ■ ■ ■ , H}. If the atom we consider 
is with a longer lifetime and under weak pumping (small 
Rabi frequency) such as the usual case in optical fre- 
quency, the time scale of the atomic evolution will be 
always longer than the time scale of the memory func- 
tions in this case. Therefore under such assumptions it 
should be legitimate to simply apply the zero-th order 
perturbation terms (zeroth-order Born approximation). 
This is equivalent to use the equal time operator prod- 
ucts to replace the two-time operator products. With 
these approximations and by using the identities of Pauli 
matrices, the generalized optical Bloch equations can be 
reduced into the following form: 



a-{t) = i—a z (t)e 



iAt 



dt'G{t-t')<7-(t') + n-{t) 



■ , s ft , s 

a+(t) = -i-a z (t) 



iAt 



(11) 



(12) 



dt'G c (t - t')a+(t') + n+(t) 
& x (t) = in{a^{ty^ - a+(t)e- lAt ) (13) 
dt'[G(t - t 1 ) + G c {t - t')](l + a z (t')) + n z (t) 



This approximation should be valid as long as the time 
scale of the memory function is still much shorter than 
the time scale of the atomic response (i.e., the inverse 
of the Rabi frequency and the decay rate). However, 
in contrast to the Markovian approximation, we do not 
approximate the memory function by a ^-function, but 
instead still keep its finite response characteristics in or- 
der to include the effects due to the non-uniform den- 
sity of states near the bandedge. Although such a sim- 
ple approximation includes only one portion of the non- 
Markovian nature of the problem, it should still be quite 
valid for the fluorescence spectrum calculation consid- 
ered in the present work since here the memory function 
time scale is typically much shorter than the atomic re- 
sponse time scale. We have also checked the validity of 
the zeroth-order approximation used here by carrying out 
the derivation in which the first order expansion term is 
also included (first-order Born approximation) |l?|. 
Theoretically the fluorescence spectrum can be calcu- 




FIG. 2: Resonance fluorescence spectrum near the bandedge 
at constant Rabi frequency. 



lated by taking the Fourier transform of the first or- 
der correlation function of the atomic dipole moment 
operator. When expressed in terms of the atomic cor- 
relation functions in the frequency domain, the expres- 
sion for the fluorescence spectrum is given by: S(lo) cx 
(a + (uj)a-{-u)) R . 

Based on the above formula, in Fig. [21 we plot the 
resonance fluorescence spectra at a constant Rabi fre- 
quency when the atomic transition frequency uj a is near 
the bandedge, uj c . All the frequency parameters used in 
the calculation are expressed in terms of the normalized 
frequency unit (3 defined above and are labeled in the fig- 
ure. It can be noted that the separation of each adjacent 
peaks is also determined by the Rabi frequency as the 
case in free space. When the atomic transition frequency 
is far away from the bandedge {oj a 3> u> c ), the normal 
resonance fluorescence spectrum of Mollow's triplets is 
obtained just as expected. However, as the atomic transi- 
tion frequency moving toward the bandedge, the profiles 
of incoherent scattering processes become sharper and 
sharper as there are fewer and fewer DOS available. The 
profile in the lower frequency is not only suppressed but 
also asymmetrical due to the existence of the bandgap 
and its residual profile exhibits a sharp edge as shown 
in the inset of the figure. It should also be noted that 
the peak in the higher frequency is enhanced a lot as 
can be clearly seen in the figure. Eventually the peak 
in the lower frequency will be totally suppressed when 
the atomic transition frequency is moving more toward 
the bandedge. At this time the resonance fluorescence 
spectrum now only has two peaks. More detailed results 
about the properties of resonance fluorescence intensity 
spectra will be reported elsewhere. 

To observe squeezing, we need to calculate the fluores- 
cence spectra for the quadrature field component. Theo- 
retically the quadrature field operator in the 9 phase an- 
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FIG. 3: Fluorescence quadrature spectra; Solid line: in- 
phase, Si(uj), Dashed line: out-of-phase, ^(w) with lo c — 
100/3 and Q = 0.25/3. 



gle is defined as E e {t) = e l6 E (+ \t) + e^ l9 E^(t). The 
two cases corresponding to 6 — and 9 = n/2 represent 
the in-phase and out-of-phase components of the electric 
field respectively. The spectra for the quadrature field 
can be obtained by calculating the following normally 
order variance H, |jj : 

S e (u>) =< E e {uj),Eg(-uj) > (14) 
oc -[< cr_ (w)cr_ (-uj) > e~ 2ie + < a + (w)a-(-w) > 
+ < a + (-~uj)a-(uj) > + < a + (-uj)a + (uj) > e 2l6 \ 

In the free space case, the in-phase quadrature S\(u) 
produces the central peak of the Mollow's triplet while 
the out-of-phase quadrature Sz{u) produces the two side- 
peaks. The situations are different for the case of pho- 
tonic bandgap crystals. When the transition frequency is 
near the bandedge of the photonic crystals, the bandgap 
effects modify the fluorescence intensity spectrum and 
cause asymmetric spectral profiles as we have seen in the 
Fig. (J2J. As shown in Fig. ©, for the case of pho- 



tonic bandgap crystals, the in-phase quadrature not only 
contributes to the central component but also the two 
sidebands. The out-of-phase quadrature still only con- 
tributes to the two sidebands as the case in free space. 
Moreover, we find both sidebands of the in-phase quadra- 
ture now exhibit squeezing even when f2 2 > ^T 2 . The 
higher frequency peak exhibits larger squeezing as well 
as larger fluorescence intensity, Fig. 

In conclusion, by introducing the Liouvillc operator 
expansion, we have successfully overcome some of the 
difficulties associated with the non-Markovian nature 
of the problem caused by the non-uniform distribution 
of the photon states in the photonic bandgap crystal. 
The calculated results for resonance pumping have in- 




FIG. 4: In-phase quadrature spectra near the bandedge with 
lo c = 100/3 and fi = 0.25/3. 



dicated that the resonance fluorescence spectra near a 
photonic bandgap can exhibit interesting behaviors in- 
cluding the suppression and enhancement of the Mol- 
low's triplet peaks, and the squeezing phenomena in the 
in-phase quadrature spectra. It shall be very interesting 
to see if one can actually verify these predictions experi- 
mentally. 



[1] R. Loudon, "The quantum theory of light, " 2nd ed. (Ox- 
ford University Press, New York, 1983.) 

[2] R. R. Mollow, Phys. Rev. 188, 1969 (1969). 

[3] C. Cohen- Tannoudji, in "Frontiers in Laser Spec- 
troscopy", edited by R. Balian, S. Haroche, and S. Liber- 
man (North-Holland, New York, 1977.) 

[4] H. J. Kimble, et al, Phys. Rev. Lett. 39, 691 (1977). 

[5] F. Y. Wu, et al, Phys. Rev. Lett. 38, 1077 (1977). 

[6] D. F. Walls and P. Zoller, Phys. Rev. Lett. 47, 709 (1981). 

[7] L. Mandel, Phys. Rev. Lett. 49, 136 (1982). 

[8] M. J. Collett, et al, Opt. Comm. 52, 145 (1984); 



[9] D. F. Walls and G. J. Milburn, "Quantum optics," 2nd 
ed. (Springer- Verlag, Berlin, 1984.) 
[10] Z. H. Lu, et al, Phys. Rev. Lett. 81, 3635 (1998). 
[11] S. John, Phys. Rev. Lett. 58, 2486 (1987). 
[12] E. Yablonovitch, Phys. Rev. Lett. 58, 2059 (1987). 
[13] S. John and J. Wang, Phys. Rev. Lett. 64, 2418 (1990). 
[14] S. John and T. Quang, Phys. Rev. A 50, 1764 (1994). 
[15] S. Y. Zhu, et al, Phys. Rev. Lett. 79, 205 (1997). 
[16] S. Y. Zhu, et al, Phys. Rev. Lett. 84, 2136 (2000). 
[17] M. Lewenstein, et al, Phys. Rev. Lett. 59, 775 (1987). 



